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Abstract 

In this paper, a new scheme of arbitrary high order accuracy in both space 
and time is proposed to solve hyperbolic conservative laws. Based on the idea of 
flux vector sphtting(FVS) scheme, we split all the space and time derivatives in 
the Taylor expansion of the numerical flux into two parts: one part with positive 
eigenvalues, another part with negative eigenvalues. According to a Lax-Wendroff 
procedure, all the time derivatives are then replaced by space derivatives. And 
the space derivatives is calculated by WENO reconstruction polynomial. One 
of the most important advantages of this new scheme is easy to implement .In 
addition, it should be pointed out, the procedure of calculating the space and 
time derivatives in numerical flux can be used as a building block to extend the 
current first order schemes to very high order accuracy in both space and time. 
Numerous numerical tests for linear and nonlinear hyperbolic conservative laws 
demonstrate that new scheme is robust and can be high order accuracy in both 
space and time. 

Keywords. Flux vector splitting scheme, arbitrary high order accuracy , space and 
time, hyperbolic conservative laws 

1 Introduction 

In the recent years, numerical schemes of high order accuracy for hyperbolic conservative 
laws have attracted much attention. In the 1980s, thanks to the efforts of Harten [7], 
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Van Leer [18] and Roe [I2] et al. (only a partial list is given here), schemes of second 
order accuracy tend to maturity. Later, ENO/WENO schemes show a possibility to 
construct the schemes of arbitrary high order accuracy in space mn- And since the 
end of the 1990s, designing schemes of high order accuracy (greater than second order 
accuracy) becomes one of the most important issues in CFD. In practical simulations, 
ENO/WENO schemes usually couple with the multi-stage Runge-Kutta method for 
time evolution. In implementation , the third order TVD Runge-Kutta method is often 
used [1]. This is because higher than third order Runge-Kutta methods will become 
complicate. For example, fourth order TVD Runge-Kutta method needs to save all the 
variables in intermediate stages, which will definitely introduce more computational 
complexity. What’s worse, higher than fourth order Runge-Kutta method will meet 
the so-called Butcher barrier [3], which the number of stages will greater than the order 
of accuracy (e.g., fifth order Runge-Kutta method need six stages). But third order 
Runge-Kutta method will also lead to accuracy barrier, for example, the fifth order 
WENO scheme couples with the third order TVD Runge-Kutta method can only achieve 
third order accuracy [IS]. To obtain ideal convergence rate, one has to reduce the CEL 
condition number, which will increase the computational costs. 

At the first decade of this century, many researchers made great efforts to avoid the 
above drawback. The first milestone may belong to the ADER (Arbitrary DERivative 
in space and time) scheme, which is a Godunov approach based scheme of arbitrary high 
order accuracy in both space and time [15]. The idea of ADER scheme can go back to 
the GRP scheme [2] , which solves the generalized Riemann problem with piecewise lin¬ 
ear initial values instead of the conventional Riemann problem with piecewise constant 
initial values. The GRP scheme is sophisticated, but may become quite complicated in 
the case of higher than third order accuracy. By introducing a linearization technique, 
the ADER scheme transfers a generalized Riemann problem into a conventional Rie¬ 
mann problem coupled with a series of linear hyperbolic problems. This transfer process 
much simplihes the computational process of the original GRP scheme and makes it 
possible to construct schemes of arbitrary high order accuracy both in space and time. 

The ADER scheme can be divided into two camps, the one is state-expansion na, 
another one is flux-expansion [HI- The former can not be integrated explicitly in time, 
so a Gaussian quadrature is necessary for high order accuracy. On the other hand, 
the latter can be integrated exactly in time, the time derivatives of the numerical flux, 
however, have to be given explicitly, and it could be very tedious when the formulation 
of the numerical flux is complex. If one chooses a nonlinear Riemann problem solver 
such as HLLG, it is not trivial to deduce the formulation of high order derivative terms. 
That is why the flux expansion version ADER is more difficult to implement than the 
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state expansion one. 

Besides ADER scheme, HGKS (high-order accurate gas-kinetic scheme) is another 
important scheme of high order accuracy pun], which is a gas kinetic scheme. HGKS 
is based on the Boltzmann equation, then the time derivatives in the Taylor expansion 
of the particle distribution function can be replaced by space derivatives directly. Gom- 
paring to the ADER scheme, the numerical flux in HGKS can be obtained without a 
linearization process. HGKS can be also integrated exactly in time and made to be of 
arbitrary high order accuracy both in space and time. We should point out, however, 
that the formulation of HGKS may also become quite complex in the high order accu¬ 
racy case. Although the recursive techniques simplify the procedure of deduction in the 
high order HGKS, the computational costs can not be reduced. In fact, the computa¬ 
tional costs in GPU time for HGKS may greater than those for WENO with Godunov 
Riemann problem solvers, see HU. Another drawback for HGKS is that there are some 
spurious velocity and pressure oscillations near the contact discontinuous when HGKS 
is used to solve the Euler equations. This phenomenon is found and carefully analyzed 
by the authors in [4]. A remedy for the GKS coupled with the Runge-Kutta method 
is also provided in their articles HE]. However, how to construct an oscillation-free 
HGKS (greater than second order accuracy) is still unsolved. 

In recent hfteen years, the ADER scheme and HGKS have provided a new approach 
to construct schemes of high order accuracy. In this article, we present a new simple 
version to construct a high order scheme and propose thus a numerical scheme of ar¬ 
bitrary high order accuracy in both space and time. Different from the ADER scheme 
and HGKS, our new scheme is based on the flux vector splitting method, and is much 
easier to implement and robust. We call the new scheme HE VS (arbitrary high order 
flux vector splitting scheme). We shall give a number of numerical tests to demonstrate 
the efficiency and accuracy of HE VS. 

This paper is organized as follows. In Section 2, we present the framework of 
numerical schemes of arbitrary high order accuracy in both space and time, while in 
Section 3 we give the construction of our new scheme HFVS in details. In Section 4 
we test HFVS in accuracy by a number of numerical examples of linear and nonlinear 
conservative laws in ID and 2D, and compare the numerical results computed by HFVS 
and those by WENO in both accuracy and computational costs. Finally, we draw the 
conclusions. 
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2 Framework of numerical schemes of high order 
accuracy in both space and time 


As mentioned in the previous section, the ADER scheme is a Godunov scheme based 
approach, while HGKS is a Boltzmann equation based approach. In our opinion, the 
flux-expansion version of the ADER scheme and HGKS can be classified in a unified 
framework of high order schemes in both space and time. In fact, Toro’s original 
article [1^ implies the possibility. Here we want to make the framework more clearly. 

Let us start with a one-dimensional hyperbolic system of conservation laws in the 
form: 


dW dF(W) , , 

To simplify the presentation, we consider here the ID Euler equations of compressible 
flows in this section, although both the ADER scheme and HGKS are constructed for 
general hyperbolic conservation laws. For the the ID Euler equations, W and F in fl2.ip 
are given by 


W={p,pU,pEf, 

F{W) = {pU, pU^ + P, pEU + PUY, 


( 2 . 2 ) 


where, p, P, P, P are the density, fluid velocity, specific total energy and pressure, 
respectively. The system can be closed by adding the equations of state (perfect gases): 


P = (7 - l)p(£ - 

Integrating (12.H over a space-time computational cell [x•_!,x-^i] x we 


have 


r 1 

1 

/,» j 

A-1 ' 


dt 


+ 


dx 


dxdt = 0. 


(2.3) 


Denote Ax = [Xj_i, x^._,_i] and At = [P,t"'+^], we approximate fl2.3p by the classical 
finite volume method (FVM) 


n 


n+l ^ n 

W, = IT. - — 
^ Ax 


7T1+1 




't" 



(2.4) 


where Wj is the cell average of W(x,F) (x G [xj_i,Xj^i]) and Pj+i is the numeri¬ 
cal flux. To make the scheme fl2.4p high order accuracy, the key point is to suitably 
construct the numerical flux which should be in high order accuracy. 

The procedure to construct a numerical scheme with high order accuracy in both 
space and time can be divided into three steps. 
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Step I. Construct F (O’*") by the Taylor expansion of the numerical flux as follows. 

N rsU h 

F,,.(r) = F,,.(0P + 5 ; (2.5) 

k=l 

Step II. Replace the time derivatives by spatial derivatives. 

Step III. Calculate the space derivatives by certain reconstruction techniques. 

In the above procedure, a WENO reconstruction technique is used in Step III for 
the ADER scheme and HGKS, while very different techniques are employed in Steps I 
and II, and described in the following. 


2.1 The ADER scheme 


In the ADER scheme, an approximate Riemann solver such as HLLC is used in Step 
II; 





F. + A,- IT, - W. 


- IE, 


if 0 < Sj, 
if Sj<0< S% 
tf S*<0< 5,+i, 




if 0 > Sj, 


where 




Pk 


(Sk-Uk\ 

\Sk-U*) 


Ek + 


1 


Uk) 



Pk. 


Pki^k Pk) 


,k 


j or j + 1 


The wave speed of Sj, Sj + 1 and S* can have many choices, see [T6] . 

Then, the time derivatives of numerical flux in (12.bh can be given analytically, for 
example, as follows. 


dF 

dt 

d'^F 

W 


(pU),, ipU), U + (pU) Ut + Ft, Ut ipE + P) + U{pE + P),J , 
\pU)u, {pU)uU + 2{pU)tUt + {pU)Uu + Pu, UuipE + P) 


T 


+2Ut{pE + P)t + U{pE + P)u 


According to the construction procedure of the Lax-Wendroff scheme, the time 
derivatives in the above formulations can be replace by the spatial derivatives. 

Remark 1. The procedure can be carried out for any given order. However, it 
is easy to see that these formulations may become very complicate in the case of very 
high order accuracy , which will increase the complexity in implementation. 


5 















2.2 HGKS 


Different from the ADER scheme, HGKS is based on the Boltzmann-type eqnation, 
snch as the BGK model: 

r . r 9 - f 

Jt + uU = -, 

T 

where / is the particle velocity distribntion fnnction and g is corresponding eqnilibrinm 
distribntion. Both / and g are fnnctions of space x, time t, the particle velocity m, and 
the internal variables ^ = (^ 1 ,^ 2 , • ■ • which has K degrees of freedom. 

The relationship between the nnmerical flnx and the particle distribntion fnnction 
is 

i (r) = (r)d5. 

Here T is the moment vector 

= (1) 2 ^ 

and the internal variable is eqnal to .^2 + ''' + ilk- 

Hence, the time derivatives of nnmerical flnx can be calcnlated by 


In view of the BGK model, the time derivatives can be evalnated directly, for ex¬ 
ample. 


/,+i(0+) = (1 - r(a'n + A^))H{u)g^ + (1 - r(a^n + A^))(l - H{u))g^, 

^/j+i(0+) = ^go - T{a^ + b)u^go - t(A‘^ + B')go - 2 t(C + aA)ugo 
+ [~ + A^ + r(G + a^A^)u + r((a')' + b^)u^]H{u)g^ 

+ [-- + A'- + t{C^ + a^A^)u + + b^y]{l - H{u))g\ 

T 

11 1 

77;^/,+i(0^) = —^90 - o.ugo H —Ago -|- (a^ -|- b)v?gQ + (A^ -|- B')gQ + 2{C + aA)ugo 

oB “' 2 T‘^ T T 

+ [^ - \alu + A') + -a^u - 2{& + a}A^)u - {{a}f + b^)u^]H{u)g^ 

T T 

+ [lr - haku + AG + -aku - 2(G" + akAAu - ((aG^ + _ Hiu))gk 

T T 

Here go is the Maxwellian distribntion at the cell interface g^ and g^ are the left and 
right limits of the Maxwellian distribntion at the cell interface, repsectively. a,b,A,--- 
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are related to the space and time derivatives of the the Maxwellian distribution, see 
[To] for the detailed dehnition. 

Remark 2. Theoretically, HGKS can achieve arbitrary high order accuracy. How¬ 
ever, due to its complexity in the expression, only third and fourth order schemes are 
often used in simulations p. fTO] . 

In the following section, we want to present a simpler high order scheme, which is 
different from the ADER scheme and HGKS. 


3 HFVS: Arbitrary high order flux vector splitting 
scheme 

From the last section, we can see that the complexity of the flux expansion version 
of the ADER scheme and HGKS actually comes from the procedure of evaluating the 
time derivatives of the numerical fluxes. So, if one can End a numerical flux, the time 
derivatives of which can be calculated very easily, then the complexity of the associated 
numerical scheme can be largely reduced. This motivates us to construct our scheme. 

As is well-known, in comparison with the traditional Godunov scheme, flux vector 
splitting (FVS) schemes are simpler in construction and coding. Thus, we try to use 
the FVS approach, instead of the Godunov approach used in the ADER scheme. Our 
proposed scheme is constructed in the same framework as in Section 2, and also consists 
of three phases: constructing spatial derivatives, leading terms and time derivatives. 


3.1 Spatial derivatives 


In order to diminish possible spurious oscillations, we should use conservative recon¬ 
struction techniques such as ENO/WENO to evaluate spatial derivatives. 

Suppose that a smooth function over a cell should be reconstructed, which is of the 
form: 

- d^W 

W{x) = + X (3.6) 

k=l 

where Wj is the cell known averaged value of W over cell j, d^W/dx^ {k = 1, - , N) 

are the unknowns to be determined, and 


(P4 



^3 


1 /X — Xj\^ 

J ’ 
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In general, we can always obtain two values on cell interfaces by the ENO/WENO 
reconstruction techniques, i.e., Thus, we immediately obtain two linear alge¬ 

braic equations: 




In this case , these formulations also implies 


1 

Ax 



W{x)dx 


W 


3 - 


(3.8) 


When the dimension = 1, the above system is uniquely solvable. 

For iV > 1, supplementary information should be added in order to have a solution 
of the system fl3.7p . A simple way of giving the supplementary information is to employ 
the cell averaged values. For example, for N = 2 we use 


1 

Ax 



W{x)dx 


Wj. 


(3.9) 


Coupling fl3.9p with (13.71) . we obtain a system of three independent linear equations for 
the three unknowns, which is uniquely solvable. 

Similarly, for = 4 we use 


1 

Ax 



W{x)dx 


Wi, 


1 = j - 1,J,J + 1, 


instead of (13.9p . Obviously, this choice is suitable for the case of arbitrary high order 
accuracy. 

In fact, the spatial derivatives can be given explicitly, see below. 

Second order. 



(3.10) 


Third order. 









Fifth order. 


- IV, - iow,_, + 

(sl) = K-W'j-i - ^OW, - + 30^^+. + 30W,_. ) 

(&) = + »'i+l - 2»';+l + 2Wi-l) 

(#) = + lOI^; + »'j+i - 6»'i+l - 


(3.12) 


3.2 Leading terms 


Following the idea of FVS, the flux in a computational cell can be decomposed into two 
parts: 




F 


i+i- 


Thus the numerical flux on a cell edge consists of two parts: 


h+i(0+)=Fho*) + f'hi(0*)- 

By a first-order Steger-Warming scheme [H], we have 

{W) = A^W, 


(3.13) 


(3.14) 


where = LA^R, L and R are the left and right eigenvectors respectively, are the 
diagonal matrices of positive resp. negative eigenvalues. On the other hand, we also 
have 

F{W) = AW, (3.15) 

where A = FAR is the Jacobian matrix, A = A"*" -|- A“ is the diagonal matrix of 
eigenvalues. 


3.3 Time derivatives 

By virtue of fl3.14p . fl2.ip and fl3.15p . omitting the subscripts, we obtain the time deriva¬ 
tives of two parts in fl3.13p : 


^F^ = (-A-^W^ .???? 

dt dt \ dx J 

Generally, we have 


Si‘ 


Qk I g 

dt 


dk 1 

dt^~^ 



-A^W 

dx 


A^{-Af 


dk ^ 

- W. 

dx^ 


(3.16) 
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Hence, we can easily get the numerical flux (in high order accuracy) as follows. 

fc=i 

where the spatial derivatives d^Wjdx^ {k = 1, - ■ ■ , N) are given in 13.10113.121 

Remark 3. From the above construction procedure, we can see the new scheme 
HFVS has two advantages. Firstly , it is obviously simpler than the original ADER 
and HGKS. The main reason is: the high order time derivatives can be expressed by 
a explicit and compact formulation of high order space derivatives, which makes the 
computational procedures clear and easy to implement. In addition, HFVS evidently 
used less computational costs than ADER and HGKS. This is also attributed to the 
simpler procedure to compute the high derivative terms. 

Remark 4. Although the Steger-Warming scheme is used as the leading term in 
the construction of our scheme, we should point out that the leading term in fl3.17p 
can be arbitrary Riemann solver, such as the HLLC and BGK solvers. This means the 
latter terms can be used as a building block to extend the current FVS scheme to high 
order accuracy in both space and time. 

4 Numerical tests 

In this section, we present the numerical experiment results for linear and nonlinear 
hyperbolic conservative laws in ID and 2D. We call our new scheme the HFVS (arbitrary 
High order Flux Vector Splitting) scheme. The abbreviations HFVS2, HFVS3 and 
HFVS5 mean the second-, third- and fifth-order HFVS schemes, respectively, owing 
to the simplicity, all numerical schemes only use conservative variables ENO/WENO 
reconstruction techniques. And we will And they work well enough in most of the cases 
below. 

4.1 Linear hyperbolic problems 

In this subsection, we want to examine the accuracy of the new schemes by solving 
some linear hyperbolic problems. 

Example 1. Accuracy test. 

We consider the following linear convection equation 






with initial data 

Wo(x) = sin(27rx). 

We take the unit interval (0,1) as the computational region and prescribe suitable 
boundary conditions (inflow boundary condition). 

Table I shows the errors and accuracy for different schemes using a same CFL 
number 0.95 at a computational time t = 1.0. We can observe that all HFVS schemes 
can reach their designed order accuracy. 


Table 1: Errors and accuracy for WENO and HFVS 


METHOD 

N 

LiERROR 

ORDER 

L 2 ERROR 

ORDER 

L^ERROR 

ORDER 

HFVS2 

20 

0.28 


0.28 


0.31 



40 

0.57E-01 

2.34 

0.65E-01 

2.15 

0.92E-01 

1.76 


80 

0.12E-01 

2.21 

0.13E-01 

2.25 

0.19E-01 

2.24 


160 

0.29E-02 

2.07 

0.32E-02 

2.08 

0.46E-02 

2.08 


320 

0.73E-03 

2.01 

0.81E-03 

2.01 

O.llE-02 

2.01 


640 

0.18E-03 

2.01 

0.20E-03 

2.01 

0.28E-03 

2.01 

HFVS3 

20 

2.32E-01 


2.54E-01 


3.35E-01 



40 

3.92E-02 

2.56 

4.38E-02 

2.53 

6.14E-02 

2.44 


80 

5.26E-03 

2.89 

5.86E-03 

2.90 

8.27E-03 

2.89 


160 

6.62E-04 

2.99 

7.36E-04 

2.99 

1.04E-03 

2.99 


320 

8.13E-05 

3.02 

9.03E-05 

3.02 

1.28E-04 

3.02 


640 

1.02E-05 

2.99 

1.13E-05 

2.99 

1.60E-05 

3.00 

HFVS5 

20 

1.02E-01 


1.17E-01 


1.59E-01 



40 

4.75E-03 

4.42 

5.35E-03 

4.45 

7.51E-03 

4.40 


80 

1.62E-04 

4.87 

1.81E-04 

4.88 

2.55E-04 

4.88 


160 

5.14E-06 

4.97 

5.71E-06 

4.98 

8.08E-06 

4.97 


320 

1.58E-07 

5.02 

1.75E-07 

5.02 

2.48E-07 

5.02 


640 

4.93E-09 

5.00 

5.48E-09 

4.99 

7.75E-09 

5.00 


4.2 Nonlinear hyperbolic problem 

In this section, we continue to present the numerical results of ID and 2D Euler equa¬ 
tions of compressible flows. 
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4.2.1 ID Euler equations 


The Governing eqnations of ID Euler equations can see 12.21 In this section, 7 = 1.4 
,and all numerical schemes use an optimal CFL number 0.95. 

Example 2. Shu-Osher shock turbulence interaction problem |13| . 
the computational domain is [— 1 , 1 ], and the initial condition is 



(3.857143, 2.629369,10.333333), a; < -0.8 
(1 + 0.2 sin 57 rx, 0.0,1.0) ,a:>—0.8 


Figure 1: Example 2. density prohles for different schemes. 


FigsU] and Figsj2] shows the density , pressure and internal prohles computed by 
different schemes t = 0.47 with mesh of 200 cells. And a reference solution is computed 
by the third order WENO scheme on a hue mesh with 2000 cells. We can observe that 
all schemes can resolve the complex solution of the governing equation. In addition, 
HFVS2 introduces over diffusion, and HFVS3 improves the resolution, while HFVS5 
gives a sharp prohle. This also demonstrates the advantages of high order scheme. 
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Figure 2: Example 2. internal energy profiles for different schemes. 

Example 3 Woodward-Colella blast wave problem| ll9] . The computational 
domain is [0,1] with reflected boundary conditions on both sides, and the initial condi¬ 
tion is 


{ (1,0,1000), 0 < X < 0.1 

(1,0, 0.01), 0.1 < X < 0.9 
(1,0,100), 0.9 < X < 1.0 

This is problem with high ratio of pressure and the solution includes strong shock 
waves, contact discontinuities, rarefaction waves and their interactions. So it is a very 
challenging problem for high order schemes. 

Figl3]and FigJU shows the numerical results at computational time t = 0.038 with 
mesh of 800 cells. The reference solution is computed by third order WENO on 10000 
cells. All HFVS schemes can resolve all the waves of the flow. As the previous test prob¬ 
lem, HFVS2 gives poor solution, HFVS3 greatly improve the resolution, and HFVS5 
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can give a sharp picture. The numerical results can demonstrate the robustness and 
accuracy of HFVS. 


Figure 3: Example 3. density prohles for different schemes. 


4.2.2 2D nonlinear problems 

In this section, we present the numerical results of 2D Euler equations: 

dW dF{W) dG{W) 

dt ^ dx dx 


where 

W = {p,pU,pV,pEf 


F{W) = {pU, pU^ + P, pUV, pEU + PUf 
G{W) = {pV, pUV, pV^ + P, pEV + PVf 
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Figure 4: Example 3. internal energy profiles for different schemes. 


and the equations of state is 

In the following case , 7 = 1.4. 

In this section, all numerical schemes use an optimal CFL number of 0.45. 
Example 4. 2D Riemann problem. We want to show two cases which are 
simple extensions of ID Riemann problems. 

The hrst case sets computational domain as [0,1] x [0,1]. And the boundary condi¬ 
tions are all outflow condition. The initial condition is given as 


{p,U,V,P) 


' (1.0, 0.0, 0.0,1.0) , X < 0.5 and y < 0.5 
(0.1, 0.0, 0.0, 0.1) , X < 0.5 and y > 0.5 
(1.0, 0.0, 0.0,1.0), X > 0.5 and y > 0.5 
^ (0.1, 0.0, 0.0, 0.1), X > 0.5 and y < 0.5 
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We continue to compare the numerical results of HFVS with different order accu¬ 
racy. The numerical results using the computational cells of 100 x 100 at f = 0.15 are 
presented. FigJT] shows the density contours between the range of 0.1 and 1.0. We can 
draw the same conclusions as ID case, i.e., all schemes can work well and HFVS can 
improve the accuracy and resolution if we increase the order of accuracy. 


Figure 5: Example 4. density contours for different schemes,from left to right: HFVS2, 
HFVS3, HFVS5. 50 contours were fit between a range of 0.1 to 1.0. 

The second case is a long time test. The computational domain is [0,1] x [0,1]. And 
the boundary conditions are all reflective conditions. The initial condition is given as 



(1, 0, 0,1), - 0.5) V (1/- 0.5f < 0.3 

(0.125,0,0,0.1) ,e/se 


The numerical results using the computational cells of 100 x 100 at t = 1.0 are 


presented. FigJT] shows the density contours between the range of 0.1 and 0.6. We 
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can draw the same conclusions as previous one, i.e., second order scheme gives a over¬ 
diffusion picture, third order scheme improve the accuracy, and fifth order scheme can 
give a sharp profile. It also demonstrates the efficiency of high order schemes. 


Figure 6: Example 4. density contours for different schemes,from left to right: HFVS2, 
HFVS3, HFVS5. 50 contours were fit between a range of 0.1 to 1.0. 


4.3 Further discussions 

4.3.1 The effective of different leading terms 

As mentioned before, it seems possible to use different leading terms to couple with 
the same high derivative terms in HFVS. We compare three problem solver: Steger- 
Warming , HLLC and BGK, which represent FVS, Godunov and GKS approaches. 
FigJT] shows second order results and FigJH] shows fifth order results of Shu-Osher prob¬ 
lem. We can observe that all schemes using different leading terms can work well. And 
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there are small differences among these scheme in low order case , while these differences 
almost vanish in high order case. 

These numerical results do really demonstrate the possibility: the high order deriva¬ 
tives terms of HFVS can be used as a building block to couple with any Riemann 
problem solvers. 


Figure 7: Shu-Osher problem, density profiles for HFVS2 using different leading 
terms.left : global picture, right: local enlargement. 


Figure 8: Shu-Osher problem, density profiles for HFVS5 using different leading terms, 
left : global picture, right: local enlargement. 


4.3.2 Comparison between WENO and HFVS 

In fact, the detailed comparison between WENO coupling with Runge- Kutta method 
with ADER can be found in na. It clearly show the advantages of ADER scheme. 
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Similarly, to demonstrate the efficiency of HFVS, we present the numerical results 
computed by WENO schemes coupling with Runge-Kutta methods. We should point 
out that we use the same reconstruction step in both WENO and HFVS. 

Here, WEN03 and WEN05 represent third order and fifth order WENO schemes , 
and RK3 is third order Runge-Kutta method. 

Table 2 shows the errors and accuracy of WEN03-I-RK3 and WEN05-I-RK3 using 
the same computational condition as Table 1. It is clearly, although WEN05-I-RK3 can 
get better accuracy than WENo3-|-RK3 , it can also achieve third order accuracy. 

Figj9] show the numerical results of Shu-Osher problems, which are computed by 
WENO and HFVS. The pictures show HFVS get better accuracy than WENO, no 
matter in case of third order or hfth order. This conclusion is similar to ADER and 
HGKS. 

Table 3 show the computational cost of WENO and HFVS. We can also observe 
that HFVS almost use half CPU time of WENO. The reason is also clear: In every 
computational time step, HFVS use only one reconstruction step, while WENO should 
use three reconstruction steps. Although HFVS should calculate the high order terms, 
it is cheaper than reconstruction step. 

Table 2: Errors and accuracy for WENO 


METHOD 

N 

LiERROR 

ORDER 

L 2 ERROR 

ORDER 

LqoERROR 

ORDER 

WEN03+RK3 

20 

7.07E-01 


7.07E-01 


9.07E-01 



40 

5.25E-01 

0.429 

5.69E-01 

0.313 

7.54E-01 

0.266 


80 

1.26E-01 

2.058 

1.39E-01 

2.033 

1.94E-01 

1.958 


160 

1.77E-02 

2.831 

1.96E-02 

2.826 

2.76E-02 

2.813 


320 

2.24E-03 

2.982 

2.49E-03 

2.976 

3.52E-03 

2.971 


640 

2.81E-04 

2.994 

3.13E-04 

2.991 

4.42E-04 

2.993 

WEN05+RK3 

20 

7.04E-01 


7.04E-01 


7.07E-01 



40 

3.00E-01 

1.230 

3.40E-01 

1.050 

4.78E-01 

0.564 


80 

4.42E-02 

2.762 

4.92E-02 

2.788 

6.94E-02 

2.784 


160 

5.47E-03 

3.014 

6.07E-03 

3.018 

8.57E-03 

3.017 


320 

6.78E-04 

3.012 

7.53E-04 

3.010 

1.06E-03 

3.015 


640 

8.45E-05 

3.004 

9.39E-05 

3.003 

1.33E-04 

2.994 
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Table 3: CPU time for WENO and HFVS 


Method 

Steps 

Cpu time 

Ratio 

WEN03+RK3 

1341 

0.42 

1 

WEN05+RK3 

1341 

0.8268 

1.96 

HFVS3 

1341 

0.218 

0.519 

HFVS5 

1341 

0.374 

0.89 


5 Conclusions 

In this paper, based on the idea of flux vector splitting scheme , we proposed a new 
scheme of arbitrary high order accuracy in both space and time to solve linear and 
nonlinear hyperbolic conservative laws. We start from splitting all the space and time 
derivatives in the Taylor expansion of the numerical flux into two parts: one part with 
positive eigenvalues, another part with negative eigenvalues. Then all the time deriva¬ 
tives can be replaced by space derivatives, according to a Lax-Wendroff procedure. A 
state-of-art WENO reconstruction polynomial is used to calculate the space derivatives 
. The new scheme is easy to implement, which will be very attractive for large CFD 
software. 

In addition, there is an interesting result: the procedure of calculating high order 
terms of numerical flux can be used as a building block to extend the current first order 
schemes to very high order accuracy in both space and time. 

Numerous numerical tests for linear and nonlinear hyperbolic conservative laws are 
presented to demonstrate that new scheme is robust and can be high order accuracy in 
both space and time. 
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